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Abstract: An explicit marching on-in-time scheme for analyzing transient 
electromagnetic wave interactions on ferromagnetic scatterers is described. The proposed 
method solves a coupled system of time domain magnetic field volume integral and 
Landau-Lifshitz-Gilbert (LLG) equations. The unknown fluxes and fields are discretized 
using full and half Schaubert-Wilton-Glisson functions in space and bandlimited 
temporal interpolation functions in time. The coupled system is cast in the form of an 
ordinary differential equation and integrated in time using a PE(CE)” type linear 
multistep method to obtain the unknown expansion coefficients. Numerical results 
demonstrating the stability and accuracy of the proposed scheme are presented. 
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" introduction 
¢ Applications 
° Existing methods 
¢ Motivation 


= Formulation 

e Coupled LLG and VIE 

¢ Discretization 

e Coupled matrix system 
= Numerical Results 

¢ Verification and convergence 
= Conclusion 
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TEML Applications 


=" Isolator!" 
¢ Microwave signal can only be 
transmitted in one direction 
e Protection against reverse power 
flow 
= Çirculator[2] 
e Waves can only be transmitted in 
rotation a 
¢ Switch in between transmitter AS 
and receiver 
Phase shifter 
e DC magnetic field to change the 
phase of the wave 
= Phased arrays 
e Reconfigurable Antennas 


¢ DC magnetic field for dynamic = “nttp:/www.tp 
tuning [2] http://file.sci 







FERRITE 
CYLINDERS 






NPERMANENT DIELECTRIC 
MAGNET SUBSTRATE 








Gc z, 7 
: om/neets/book11/44q.htm 
-Org/Html/3-9801086_2766.htm 
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CEML Existing Methods s 
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= Full wave simulation tools are used to design 
and characterize these systems 

= LLG should be considered while solving 
Maxwell equations for fields/currents. 

= Differential equation methods 


e LLG Is converted into a tensor and introduced as a 
constitutive relation between magnetic flux and 
field 


e LLG is introduced as an auxiliary equation 
= Integral equation methods 


e LLG is converted into a tensor and introduced as a 
constitutive relation between magnetic flux and 
field 
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SEMEL Existing Methods ‘eo? 


= Tensor method 
e Reduced number of unknowns 
e Simplified via small signal approximation 
e Accurate only around saturation point 
e Dispersive 
¢ Difficult and costly to implement in time domain 


" Auxiliary equation method 
e Flexible: Non-linearity can be modeled 
e Higher number of unknowns 
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CEMEL Motivation 
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= Differential equation methods 
e Volumetric discretization of the whole computation 
domain 
e Approximate absorbing boundary conditions (PML, 
ABC) 
e Time step size should satisfy Courant condition 


n integral Equation methods 
Discretization of scatterer 
e implicitly satisfy radiation condition 
e Time step size should resolve only the incident 
field’s maximum frequency 
e LLG is introduced as an auxiliary equation 
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Formulation (x 


= Scattering problem 


M, (r,t) : Equivalent current density 
M (r,t)=0, B(r,t)- u, 0, H(r,t) 
H'"° (r,t): Incident magnetic field intensity 
H” (r,t): Scattered magnetic field intensity 
H(r,t) : Total magnetic field intensity 
H(r, £) = H™ (r,t)+ H” (r,t) 
B(r,t) : Total magnetic flux density 
E, : Free space permittivity 


Lo : Free space permeability A 


L H™ (r,t) 
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H™ (r,t) 


{E Lho? 
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CEMEL Formulation (o 
= TD-MFVIE 
0H (r,t)=0 H(r,t)+ C{H(r,t)}— Liu, Bir,t)} 
rev 
LIX) = = ole PXE), ae 
R=|r-r' 
aval” I XE tare d 
Co = 1/ Ey Mo 
= LLG 
0,M (r,t) = an (r,4)x H (r,t) rev 
EP e © y —“—_ M_ (r,t) x | M,(r,¢)x H,(r,0)} 


VY . Damping constant 
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CEMEC Discretization 

= Expansion of unknowns 
Hr) => {LO}, £2) 
Birt) => 0 {1}, f(r) 





Half and full SWG functions : f'(r),f, (r) 
Time dependent unknown coefficients : L(z), I(t) 


A 
f" (r)=— 
„ (r) 3 





+ + 
(r-r) rey, 


+ 
k 





BID. H. Schaubert, et al., IEEE Trans. Antennas Propag., 32(1), pp. 77—85, 1984. 
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CEML Discretization 
" Spatially testing with f"(r) 
G" 0 L(t)=V™" (O+ VELA) + VEO) 


= Gram matrix 


! 


A 
Y 
wy 


( 


(Gna = f, m0) f (r)dr 


= Tested incident field 


(V"™(O},, =|) fh(r)-0,H™ (r, Oder 
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CEML Discretization 
" Spatially testing with f"(r) 
G” A L) = VIE (O+ VELO) + VEO) 


( 


= Tested scattered field 


Vi CLO}, = Def, aO LLO £0) fae 


(VECI hn = -2a f, 0L {u0 AOE, OO) fae 
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CEML Discretization (so 
= Temporal sampling with time step size Af 

L, = L(iAt) 

IL =I(iAt) 


= To accurately determine the retarded-time integrals in 
operator C{-}, L,and I should be interpolated 


L()= >) TOL, = > 7 (t-iADL, 
I(t) = LTO — KC ~ FAI), 


T (t) : Temporal interpolation functions 
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CEML Discretization 


= MOT matrix equation 
G"0,L,=Z)L,-p,Z, 1b 


Ci 


+V 4 $ ZL -£ ula y 


j-i i j—i i 


vy ixe ed 
J 
= Matrix elements 


LP Sun = J, fn) LEE OTOH, ae 
ZT yan =], OLE OT OY, de 
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CEMEL Discretization 


"= LLG — 
OQ * Gyro-magnetic ratio 
0, M, (r,t) — r (r, t) x H Ar, t) y. . Damping constant 


M, (r,t)x {M,(r,t)x H,(r,t)} 
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0) 
wy 


( 





Hae 
M (r,t)= u, B (r,t)—- H, (r,t) 





" Lossless: 0,M,(r,t)=-yM,(r,t)x H (r,t) 


= With external DC bias field 
H,(r,t)=H,,.(r)+ H(r,2) 
M (r,¢)=M_(r)+ 4, B(r,t)— H(r,7) 
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CEML Discretization 


= Under small signal approximation 


( 


O,M(r,t)=7H,,(r) x ‘B(r, 7) 


-y |M (r)+ H, (r)|x H(r,0) 
" Auxiliary equation 


Ly 0, B(r,t)—- 0,H(r,7) — yH, (r) x Ly B(r,t) 
-y [M, (r)+ H, (r)| x H(r, 2) 
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CEML Discretization 


= Spatially tested with f(r) and sampled in time 





H, G” 8,1, -G* 3, L, = H cth 1,- M,L, 


j 
= Matrix elements 


(Hpclnn = | E0) -yH 0) xf? T)dr 


V 


m 


(Mshua = | £,00)-7[M, (0) +H, (r)] xf) (nde 


m 


(Gna = f EOE dr Gn, = f, E-E r)dr 
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CEML Discretization 


= MOT matrix equation 


hh hh -l rz hf 
GƏ L = ZPL -u'zi 





j—l j—l 
inc hh —] -rz hf 
+V, + > Z,- > ly ZA, 
= [= 


vines 





= Constitutive equation 


H, G” 0,1, -G" 0, L, =H c4 1, -M,L, 
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CEMEL MOT Formulation < 
= Coupled matrix system 
Gh 4 AL Z) Lu, Z,1,+V)" 


a ab —Igvil _] 
G 4G o I, H pcth 1,- M,L, 


"= Relates time derivative of samples to samples 


= Integrated in time using predictor-corrector methods — 
PE(CE)™ 
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CEMT Explicit MOT 


° Attime step J 


Jl j- 
Determine Vi" Vet ove ty ZL > ZT, 
J J J j-i i 0 jJ-ii 

| i=l i=l 
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CEML Explicit MOT 


° Attime step J 


Determine ve 


Predict IL I } 


P J 


L = > Lp}, L is +{P} 0, L a Lig 
l= 2X I jek +{P3 0,1 í ON a 
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CEMT Explicit MOT 


( 


° Attime step J 


Determine ve 


<b G" 0 o L, 


-G^ u'G* O A, 


Predict L,I} ’ 
hh —lryhf fixed 
<b | ZYL -ag Zii +V) 
Evaluate 10,L AR H, o wd — M,L, 
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Explicit MOT 





CEML 


° Attime step J 


T L’ = =S Heh La tiaL Lia 


tn {01245} i >, Lon ne Hð y Lya 
| FACS, d; "vel 
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CEML 
° Attime step J 


paren 
V 


Pee (Lb) 
Vv 





Explicit MOT 


LY =a’ +(l-a@)L” 
J J J 


vo V v—l 
I =al, + (1 — ai, 
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CEML 

° Attime step j Gm 0 o, L, 
fh — fF 

-G 4, t o, L, 


Explicit MOT 





hh -l rz hf fixed 
Vy LZ, L- 4 Zl, +V, 


| T Hct I, -M L, 


em em 
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CEML 


° Attime step J 


Explicit MOT 
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CEML 


° Attime step J 





Explicit MOT 
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CEMT Explicit MOT 





° Attime step J 
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" Excitation is a modulated Gaussian pulse 


C4 y 


G(t, =n, cos| 27 f(t-t,) le 20” 





- Center frequency : J 

_ Effective bandwidth : fow 

— Time duration : 0=3/(27f,,) 
_ Time delay ; £,=100 


- Intrinsic impedance : 1 = 4/ 4h l £o 
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= Predictor-corrector coefficients 


— Highly accurate and stable predictor and corrector 
scheme’ 


—- Number of steps > k=22 
_~ Convergence condition : 6=10 ” 


= Temporal interpolation functions T(t) 
— Non-causal band limited interpolation functions? 


— Causality retained by complex exponential 
extrapolation scheme® 


MA. Glaser and V. Rokhlin, J. of Sci. Comput., 38(3), pp. 368-399, 2009. 
SID. S. Weile, et a/., [EEE Trans. Antennas Propag., 52(1), pp. 283-295, 2004. 
61S. Bin Sayed, et al., IEEE Trans. Antennas Propag.,63(7), pp. 3098-3110, 2015. 
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= Isotropic scatterer 


"(Or t)=yG(t—r-z/c,, fo 
0, M(r,t)= 0, H(r,t) 


wae 
Se 


` 
AN 








= 150 MHz 
= 50 MHz 
At = 0.2 ns 
r= 2M E 
=19692 | o 
æ = 0.6 


, inc (r, f) 
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CEML Numerical Results 





! 


A 
0) 
wy 


Ci 


= Isotropic scatterer 
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Numerical Results CD 
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= Ferrite scatterer with DC bias 


™ (r,t) yG(t—r-z/c,, fi 
4¢ = 2.1420 Oe M, =0.7854 Gauss 


Se eS 
DISS 
Tanaman ANa 


~ 
POSS 





, = 150 

hy = OO 

At = 0.2 ns — 
, = 19692 a 
a=0.6 a 





, inc (r, f) 
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CEML Numerical Results Ge 
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= Ferrite scatterer with DC bias 


Q 
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Conclusions 
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= Summary 


e LLG equation is used in the auxiliary form to analyze 
ferrite scatterers using integral equations 


= Future work 
¢ Further verification 


¢ Implementation of LLG with no small signal 
approximation (fully nonlinear) 


¢ Application to “real-life” microwave devices 
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